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Three, interrelated biologically-relevant examples of biased random walks are presented: (1) A 
model for DNA melting, modelled as DNA unzipping, which provides a way to illustrate the role of 
the Boltzmann factor in a venue well-known to biology and pre-medical students; (2) the activity 
of helicase motor proteins in unzipping double-stranded DNA, for example, at the replication fork, 
which is an example of a Brownian ratchet; (3) force generation by actin polymerization, which 
is another Brownian ratchet, and for which the force and actin-concentration dependence of the 
velocity of actin polymerization is determined. 



I. INTRODUCTION 

In spite of the growing recognition that physics skills - 
"scholastic rigor, analytical thinking, quantitative assess- 
ment, and the analysis of complex systems" [T] - are im- 
portant for biology [2J and pre-medical [1: students, these 
students often arrive in physics classes skeptical about 
the relevance of physics to their academic and profes- 
sional goals. To engage these students, in the 2010-2011 
academic year, the Yale physics department debuted a 
new introductory physics sequence, that, in addition to 
covering the basics - kinematics, force, energy, momen- 
tum, Hooke's Law, Ohm's Law, Maxwell's equations etc. 
- also covers a number of more biologically-relevant top- 
ics, including, in particular, probability, random walks, 
and the Boltzmann factor. The point of view of the class 
is that the essential aspect of physics is that it consti- 
tutes a mathematical description of the natural world, 
irrespective of whether the topic is planetary motion or 
cellular motion. 

The enrollment in the new sequence was approximately 
100 students. The class is evenly split between sopho- 
mores and juniors with a few seniors. The majority (80%) 
are biology majors, with 80% identifying themselves as 
premedical students, and they possess considerable bio- 
logical sophistication. In many cases, they are involved 
in biomedical research at Yale or at the Yale School 
of Medicine. In many cases too, they are involved in 
medically-related volunteer work. The major time com- 
mitment required to do justice to a rigorous physics class 
has to compete with these other obligations. Therefore, 
an important aspect of our teaching strategy is to con- 
vince these students that physics is indeed relevant to 
their goals. To this end, we determined to cover a num- 
ber of biologically-relevant topics, with which the major- 
ity of the students would have some familiarity from then- 
ear lier biology and chemistry classes. 

This paper presents three such topics, that are inter- 
related and can be treated as random walks, in the hope 
that these may be useful to others. First is DNA melting 
[3], which we place in the context of polymerase chain 
reaction (PCR). This provides a way to illustrate the 
role of the Boltzmann factor in a venue well-known to 



the students. This treatment builds on earlier sections 
of the course, concerned with random walks and chem- 
ical reaction rates, which are not described here. The 
second topic is the activity of helicase motor proteins in 
unzipping double-stranded nucleic acid (DNA or RNA, 
although we will write in terms of DNA) . Our discussion 
is based on Ref. 4 Helicase activity constitutes an el- 
egant example of a Brownian ratchet and builds on the 
earlier discussion of DNA melting. Third, we present a 
discussion of force generation by actin polymerization, 
which provides the physical basis of cell motility in many 
cases, and which is another Brownian ratchet. In this 
case, based on Ref. (5j we can determine how the velocity 
of actin polymerization depends on actin concentration 
and on load. In each of these examples, biology and pre- 
medical students in an introductory physics class see that 
a physics-based approach permits a new, deeper under- 
standing of a familiar molecular-biological phenomenon. 



II. THE BOLTZMANN FACTOR 

"The laws of thermodynamics may easily be obtained 
from the principles of statistical mechanics, of which they 
are an incomplete expression." J.W. Gibbs (§]. 

Instead of introducing thermal phenomena via ther- 
modynamics and heat engines, as might occur in a tradi- 
tional introductory sequence, following the suggestion of 
Garcia et al. [7] , we chose to assert the Boltzmann factor 
as the fundamental axiom of thermal physics. Build- 
ing upon earlier sections of the course on probability 
and random walks, this approach permits us to rapidly 
progress to physics-based treatments of DNA melting, 
unzipping of double-stranded DNA at the replication 
fork by helicase motor proteins, and force-generation by 
actin-polymerization. Specifically, we assert that, for mi- 
crostates i and j of a system, the probability (pi) of real- 
izing a microstate i and the probability (pj) of realizing 
a microstate j are related via 

Pi = e -(.U-e S )/(k B T) /j\ 
Pi 

where e, is the energy of microstate i, €j is the energy 
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of microstate j, ks = 1-38 x 10 23 JK 1 is Boltzmann's 
constant, and T is the absolute temperature. 

"This fundamental law is the summit of statistical me- 
chanics, and the entire subject is either the slide-down 
from this summit, as the principle is applied to various 
cases, or the climb up to where the fundamental law is 
derived and the concepts of thermal equilibrium and tem- 
perature clarified." R. P. Feynman on the Boltzmann fac- 
tor |H]. 

To illustrate the Boltzmann factor in a simple example, 
we consider protein folding/unfolding. Protein/unfolding 
is an example of an isomerization reaction, in which one 
chemical species alternates between different molecular 
configurations. In this case, it is important to realize 
that the folded state corresponds to a single microstate, 
but that the unfolded state corresponds to g microstates. 
This is because there is just one molecular configuration 
associated with the folded state. By contrast, the un- 
folded state can be viewed MS H random walk in space, 
and therefore corresponds to g different molecular con- 
figurations, one for each different random walk. If there 
are a total of n proteins, n u of which are unfolded, and if 
there are g possible unfolded microstates, then the prob- 
ability of realizing a particular unfolded microstate (p u ) 
is equal to the probability that a protein molecule is un- 
folded multipled by the probability that an unfolded pro- 
tein is in the particular unfolded microstate of interest, 
which is one of g equally-likely microstates: 
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There is a unique folded microstate, so in terms of n 
and the number of folded proteins, nt, the probability of 
realizing the folded microstate is simply 
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Combining EQ. [TJ EQ. |2j and EQ. [3j we find 

-(e u -e f )/(k B T) 
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(3) 
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where e u is the energy of any of the unfolded states and 
ef is the energy of the folded state. 



III. DNA MELTING 

A. DNA unzipping/zipping as a chemical reaction 

Next, we examine DNA melting, according to the 
model of Ref. [3J in which DNA melting is equivalent 
to DNA unzipping. We treat DNA zipping and unzip- 
ping as a set of isomerization reactions. To this end, we 
consider a population of identical DNA strands each of 
which contains a junction between dsDNA and ssDNA. 
Fig. [T] illustrates the reactions involving the DNA strand 




FIG. 1: Chemical reaction scheme for DNA zipping and un- 
zipping. The possible reactions of a DNA strand with i zipped 
base pairs are illustrated, either undergoing isomerization to a 
DNA strand with i — 1 base pairs or isomerization to a DNA 
strand with i + 1 base pairs. 



with i paired base pairs. This is the chemical species in 
the center. The species on the left and right are DNA 
strands with i + 1 and i — 1 paired base pairs, respec- 
tively. The relevant reaction rates are a, which is the 
zipping rate, and /3, which is the unzipping rate. When 
a > j3, the DNA zips up. When a < f3, the DNA unzips. 
As suggested in Fig. [TJ rii is the mean number of DNA 
strands with i paired base pairs, etc. 

We have previously discussed in class that how the 
concentration of chemical species changes in time can be 
described by chemical rate equations. With the help of 
Fig. [T] we are thus lead to an equation for the rate of 
change of m in terms of a, /3, rii, rii-i, and rii+\. 



~dt 



= —an.i — j3rii + arii—i + Prii+\. 



(5) 



At equilibrium, at a temperature T, on-average nothing 
changes as a function of time, so drii/dt = 0. Thus, 



= —a — p + a h p 



(6) 



The factor rii—i/rii, which is the ratio of the mean num- 
ber of DNA strands with i — 1 zipped base pairs to the 
mean number with i zipped base pairs, is equal to the ra- 
tio of the probability that a particular DNA strand has 
i — 1 zipped base pairs to the probability that it has i 
zipped base pairs. Thus, this factor is given by a Boltz- 
mann factor (c/. EQ.Bl: 



i = g e ^/(^T) 



,-{e-k B Tlng)/(k B T) 



-AG/(t s T) 



(7) 



where e is the energy required to unzip one additional 
base pair (so e is positive) and g specifies that the two 
unzipped ssDNA bases have a factor g times as many 
microstates as the single dsDNA base pair they replace. 
The last equality in EQ.[7]defines the free energy required 
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to unzip one base pair: 

AG = e- k B T In g. 



(8) 



Students are familar with AG from their chemistry 
classes. Similarly, we have 



"'+1 _ e +AG/(k B T) 



(9) 



Substituting EQ. [7] and EQ.[9]mto EQ. [§ we have 

= -a - [3 + ae- AG /^ + /3 e + AG /^ T >. (10) 



It follows from EQ.flOlthat 



a -AG/(fc B T) 



(11) 



EQ. 11 informs us that the DNA unzips, i.e. f3 > a, only 
if AG < 0, i.e. only if e — ksTlng < 0. In order for this 
condition (e — kgTlng < 0) to be satisfied, it is neces- 
sary that T > e/ (ks In <?). If we define the DNA "melting 
temperature" to be Tm = e/(fcslng), we see that the 
DNA unzips for T > Tm, while it zips up for T < Tm- 
This phenomenon is an essential ingredient in DNA mul- 
tiplication by polymerase chain reaction (PCR) , which is 
well-known to the students, and for which Kary Mullis 
won the 1993 Nobel Prize in Chemistry [9]. The first 
step in PCR is to raise the temperature above Tm, so 
that each dsDNA strand unzips to become two ssDNA 
strands. When the temperature is subsequently reduced 
in the presence of oligonucleotide primers, nucleotides 
and DNA polymerase, each previously-unzipped, ssDNA 
strand templates its own conversion to dsDNA. This dou- 
bles the original number of dsDNA strands because a new 
dsDNA strand is created for each ssDNA. PCR involves 
repeating this temperature cycling process multiple (N) 
times, with the result that the initial number of dsDNA 
molecules is multiplied by a factor of 2 N . Thus, initially 
tiny quantities of dsDNA can be hugely amplified, and 
subsequently sequenced. 



B. DNA unzipping/zipping from a random-walk 
point of view 



readily calculate the mean displacement of the ssDNA- 
dsDNA junction in a time At: 



Axj = b(/3 - a) At, 



(12) 



where zipping corresponds to a negative displacement of 
the ss-to-ds junction Since the mean of the sum of n 
identically-distributed, statistically-independent random 
variables is n times the mean of one of them (which 
students learned earlier in the course), then in a time 
t = nAt the mean displacement of the ssDNA-dsDNA 
junction is 



t 



Axj = —b(/3-a)At = b(P-a)t. (13) 



The corresponding drift velocity of the ssDNA-dsDNA 
junction is 

v 3 = b{P -a) = ba{e- AG ^ kBT ^ - 1) (14) 

This is the drift velocity of a dsDNA-ssDNA junction 
in terms of the zipping up rate (a) and the unzipping 
rate (/3), or the zipping rate (a) and the unzipping free 
energy (AG). We will come back to this result below, 
but we note now that EQ. [14] is appropriate only when 
the junction is far from a helicase. 

As defined in EQ. [8] AG is the change in free energy 
that occurs when one additional base pair is unzipped. 
Thus, as far as this expression for AG is concerned, the 
final, "product" state is the unzipped state, and the ini- 
tial, "reactant" state is the zipped state. Thus, unzipping 
corresponds to the forward direction of the reaction. We 
may make contact with what students have learned in 
chemistry classes, namely that a reaction proceeds for- 
ward if AG is negative, by pointing out that EQ. 
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forms us that the unzipping reaction proceeds forwards 
(i.e. that Vj > 0) only for AG < 0, exactly as we are 
told in chemistry classes. Here, though, this result is de- 
rived from a more basic principle, namely the Boltzmann 
factor. 



IV. HELICASE DNA-UNZIPPING ACTIVITY 



A. Helicases 



It is also instructive to view DNA zipping/unzipping 
as a biased random walk, which students have previously 
studied in the class. In this context, if we consider a 
dsDNA-ssDNA junction, the probability of zipping up 
one base pair in a time At is aAt, and the probability 
of unzipping one base pair in a time At is ft At. For 
small enough At it is reasonable to assume that the only 
three possibilities are (1) to zip up one base pair or (2) to 
unzip one base pair or (3) to not do anything. Therefore, 
since probabilities sum to unity, we must have that the 
probability to do nothing is 1 — aAt — (3 At. Given these 
probabilities, and the length of a base pair, b, we may 



Helicases [TD_ are a class of motor proteins (a.k.a. 
molecular motors), which perform myriad tasks in the 
cell by catalyzing ATP-to-ADP hydrolysis and using the 
free energy released in this reaction to do work. The im- 
portance of helicases may be judged from the fact that 
4% of the yeast genome codes for some kind of helicase. 
One of their roles is to unzip dsDNA and/or dsRNA. 
Thus, helicases play an indispensible role in DNA repli- 
cation, for example. To engage students in this topic, we 
start by showing a number of online movies illustrating 
the DNA-unzipping activity of helicase motor proteins at 
the replication fork |111 H"2"| . These movies also present 
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an opportunity for active learning in which we ask stu- 
dents to discuss with their neighbors what is misleading 
about the videos. The essential point is that, wonderful 
as they are, the videos suggest that everything proceeds 
deterministically. By contrast, as we will discuss, all of 
the processes depicted are actually random walks, but 
with a drift velocity that corresponds to their progress. 
We also point out that, on the medical side, Werner syn- 
drome, which involves accelerated aging, is caused by a 
mutation in the WRN gene which codes for the helicase 
WRN [T3]. 



B. Brownian ratchet mechanism of helicase activity 




unction 



One proposed mechanism for how helicase unzips DNA 
is as follows. The helicase steps unimpeded on ssDNA to- 
wards a ss-to-ds junction, until it encounters the junction, 
which then blocks its further progress, because the heli- 
case translocates only on ssDNA. However, at the junc- 
tion, there is a non-zero probability per unit time for 
the junction to thermally unzip one base pair, because of 
the Boltzmann factor. It is then possible for the helicase 
to step into the just-unzipped position. If the helicase 
does this, the DNA is prevented from subsequently zip- 
ping back up again. In this way, the junction is unzipped 
one step. Repeating this process many times leads to 
the complete unzipping of the DNA. Because this mech- 
anism relies on random Brownian motions to both unzip 
the DNA and to move the helicase into the just-unzipped 
position, the helicase is said to be a Brownian ratchet [5] , 
analogous to Feynman's thermal ratchet |14j . 



C. Helicase translocation from a random-walk 
point of view 

Just as the motion of the ss-to-ds junction may be con- 
ceived as a random walk, so may be the translocation of 
the helicase on ssDNA. In this case, the probability of 
the helicase stepping one base pair towards the junction 
(+6) in a time At is k + At, and the probability stepping 
one base pair away from the junction (—6) in a time At is 
k- At, where fc+ and fc_ are the rate of stepping towards 
the junction and the rate of stepping away from the junc- 
tion, respectively. Since probabilities sum to unity, and 
we assume that the only three possibilities in a small time 
At are to step towards the junction one base pair or to 
step away from the junction one base pair or to not do 
anything, we must have that the probability to do noth- 
ing is 1 — k+At — k-At. Given these probabilities, and 
the length of a base pair, 6, we may calculate the mean 
displacement of the helicase in a time At: 



Ax h = b(k + - fc„)Ai. 



(15) 




FIG. 2: Schematic of helicase translocation on ssDNA and 
DNA zipping, showing the relevant parameters. In the con- 
text of this figure, the coordinate system used in our dis- 
cussion takes the ^-direction to increase towards the right, 
so that a and fc_ correspond to motion in the negative re- 
direction, and P and k+ correspond to motion in the positive 
x-direction. 



mean of one of them, then in a time t = nAt the mean 
displacement of the helicase is 

x h = nAx h = -^ t ^x h = ^b{k + -k_)At = b(k + -k_)t. 

(16) 



The corresponding drift velocity of the helicase is 
v h = b(k + -k-). 



(17) 



This is the drift velocity of a helicase in terms of 
the stepping-towards-the-junction rate (fc+) and the 
stepping-away-from-the-junction rate (fc_). Just like EQ. 
[TT] EQ. [17] is appropriate only when the helicase is far 
from the junction. 

An important additional point, concerning helicase 
translocation on ssDNA, is that, as we saw in EQ. [TTJ 
the ratio of forward and backward rates is given by a 
change in free energy. Thus, for helicase stepping we 
must expect, in analogy with EQ. [TT] that the ratio of 
stepping rates is given by 



+ = e AG'/(k B T) 



(18) 



Since the mean of the sum of n identically-distributed, 
statistically-independent random variables is n times the 



where AG' is a free energy change. But what free energy 
change? The answer can be gleaned from the observa- 
tion that helicases, and motor proteins generally, can be 
thought of as enzymes, which catalyze ATP-to-ADP hy- 
drolysis, which is coupled to the helicase's translocation. 
It follows that AG' in EQ. 18 corresponds to the free en- 
ergy difference between ATP and ADP. (Note that AG', 
as specified in EQ.[l8] must be positive, in order to ensure 
that k + > fc_ so that the helicase translocates on ssDNA 
preferentially towards the ssDNA-to-dsDNA junction.) 
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D. Clash of the titans 

So far, we have considered the situtation when the ds- 
to-ss junction and the helicase are far apart. To deter- 
mine how helicase unzips dsDNA. it is necessary to de- 
termine what happens when these two objects come into 
close proximity, given that they cannot cross each other. 
To elucidate what happens in this case, we show to the 
class a simple Mathematic Demonstration that simulates 
these two non-crossing random walks [T5]. The simula- 
tion treats both the location of the helicase and the lo- 
cation of the ds-to-ss junction as random walks. At each 
time step within the simulation, the helicase ordinarily 
steps in the positive cc-direction, towards the junction, 
with probability k + At and in the negative x-direction, 
away from the junction, with probability k_At, while 
the junction ordinarily steps in the positive x-direction, 
zipping up one step, with probability aAt and in the 
negative ^-direction, unzipping one step, with probabil- 
ity /3At. However, in the simulation, if the helicase and 
the junction are neighbors, neither one is permitted to 
step to where it would overlap with the other. Thus, 
the helicase and the junction cannot cross. An exam- 
ple of the simulational results is shown in Fig. [3] where 
the orange trace represents helicase location as a func- 
tion of time and the green trace represents the location 
of the ssDNA-to-dsDNA junction as a function of time. 
Evidently, the helicase and the junction track together, 
implying that they have the same drift velocity. For the 
parameters of this simulation, the helicase translocates in 
the same direction as it would in the absence of the junc- 
tion. By contrast, the junction's direction is opposite its 
direction without the helicase. Thus, for these param- 
eters, the helicase indeed unzips dsDNA. Using sliders 
within the Mathcmatica Demonstration, which is readily 
accessed via any web browser, students can explore for 
themselves the effects of varying a, /3, fc+, and fc_. 

To incorporate analytically the fact that the helicase 
and the junction can not cross, we introduce the proba- 
bility, P, that the helicase and the junction are not next 
to each other. The dsDNA-to-ssDNA junction can only 
zip up if the helicase and the junction are not next to 
each other. Therefore, we reason that EQ . [12] should be 
modified to read 

Axj = b{/3 - aP)At. (19) 
Similarly, EQ . fl5] should be modified to read 

Ax h =b(k+P -k_)At. (20) 
It follows that the drift velocities are modified to read 



Vj = b(j3-aP), 



and 



b(k+P~k^). 



(21) 



(22) 




4000 
Number of steps 



FIG. 3: Simulation, implemented as a Mathematica Demon- 
stration, of a helicase motor protein unzipping double- 
stranded DNA, according to the Brownian ratchet model of 
helicase activity. The orange trace represents helicase translo- 
cation on ssDNA (left) for k+ = 0.79 and fc_ = 0.07. The 
green random walk represents the position of a ssDNA-to- 
dsDNA junction for a = 7.2 and /3 = 2.45. In both cases, 
At = 0.01. The blue line is EQ. [2§] The two random walks 
start at in the case of the junction, and at -2 in the case of 
the helicase. 



However, from the simulation, it is also clear that, 
while the helicase is unzipping DNA, the helicase and 
the junction must have the same drift velocity i.e. 



or 



b{k+P-k_) =b{f3-aP). 
We can solve this equation to determine P: 



P = 



k+ + , 



(23) 



(24) 



(25) 



Furthermore, we can use this expression for P to deter- 
mine the drift velocity at which the helicase unzips the 



dsDNA by substituting into EQ. 22 Setting Vj = Vh = v, 
we find 



v = b 



fc_ 
fe+ 



(26) 



EQ. |26| represents the velocity at which helicase unzips 
dsDNA according to the Brownian ratchet mechanism. 
The numerator in EQ. [26] is the difference of two rate 
ratios. It follows, using EQ. [TT] and EQ. [18] in EQ. [26] 
that 



-AG/(k B T) _ p -AG'/(k B T) 




(27) 



EQ. [27] informs us that whether or not helicase un- 
zips dsDNA depends solely on whether AG" > AG or 
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not. For AG' > AG, the drift velocity of the helicase- 
phis-j unction is positive, corresponding to the helicase 
unzipping the dsDNA. In fact, for one base pair, we 
have AG ~ 3fcsT, while for the hydrolysis of one ATP 
molecule, we have AG' ~ 16fcgT, so indeed the helicase 
has plenty of free energy to do its work. In fact, energet- 
ically, one ATP hydrolysis cycle could unzip up to about 
5 base pairs. 

In fact, beautiful, single-helicase experiments jTS] sug- 
gest that the simple Brownian ratchet mechanism of he- 
licase DNA-unzipping activity, presented here, should be 
refined by incorporating both a softer repulsive potential 
between the helicase and the ds-to-ss junction than the 
hard-wall potential implicit in our discussion, and suit- 
able free energy barriers between different microstates 
of the helicase and junction [3J [TB]. Appropriate 
choices of the potential and the barriers permit the heli- 
case to unzip dsDNA faster than would occur in the case 
of a hard-wall potential. A force (/), that tends to un- 
zip the DNA, can be incorporated by replacing a with 
ae -fb/(k B T) and p with p e fb/(k B T)_ 

V. FORCE GENERATION BY ACTIN 
POLYMERIZATION 

A. Actin polymerization is a Brownian ratchet 




The mechanism by which actin or tubulin polymeriza- 
tion exerts a force also a Brownian ratchet [T7J [TS] 
and is schematically illustrated in Fig. [4] In the case of 
a load, /, applied to a cell membrane, the cell membrane 
is in turn pushed against the tip of an actin filament, 
which usually prevents the addition of an additional actin 
monomer (G-actin) of length a to the tip of the actin fil- 
ament (F-actin). However, with probability specified by 
a Boltzmann factor, the membrane's position relative to 
the tip, z, occasionally fluctuates far enough away from 
the filament tip (z > a) to allow a monomer to fit into 
the gap. If a monomer does indeed insert and add to the 
end of the filament, the result is that the filament and 
therefore the membrane move one step forward, doing 
work against the load force. Repeating this many times 
for many such filaments gives rise to cell motility against 
viscous forces. In class we also show movies showing cells 
moving as a result of actin polymerization |19) . and Lis- 
teria monocytogenes actin "rockets" |20| . 

B. Actin polymerization as a biased random walk 

Similarly to EQ. [T2|and EQ. [15] we can write down an 
expression for the mean displacement of the filament tip 
in a time At in the absence of a nearby membrane: 

Ax = a(ck + - fc_)Ai (28) 

where c is the concentration of G-actin, a is the length 
of an actin monomer, k + is the actin on-rate, and fc_ is 



FIG. 4: Cartoon illustrating how actin polymerization can 
do work against a load (/), applied to a membrane against 
which the polymerizing actin filament (F-actin) abuts. Only 
if the gap, z, between the tip of the actin filament and the 
membrane exceeds the length, a, of a G-actin monomer is it 
possible for the filament to grow. 

the actin off-rate. However, if the membrane is nearby, it 
is only possible to add an actin monomer if the distance 
between the filament and the membrane is greater than a. 
Assuming that the time-scale for membrane fluctuations 
is much faster than that for adding actin monomers, if 
the probability, that the membrane-filament tip distance 
is greater than a. is P, then EQ.[32]is modified to read 

Ax = a(ck + P - fc_)At, (29) 

and the drift velocity of the tip is 

v = a(ck + P - k-). (30) 

But application of the Boltzmann factor informs us that, 
when the force on the membrane is /, the probability 
that the gap is greater than a is 

P = e -/./M ) (3i) 

so that 

v = a{ck+e- fa/{kBT) - k-). (32) 

This is the force-velocity relationship for an actin |TB] or 
tubulin [T7] filament. Although the load, /, is applied 
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to the membrane, the velocity is constant. Therefore, 
according to Newton's third Law, as the students know, 
there can be no net force on the membrane. We may de- 
duce that the load is balanced by an equal and oppposite 
force, generated by the polymerization ratchet. 

VI. CONCLUSIONS 

Three, interrelated biologically-relevant examples of 
biased random walks were presented. First, we presented 
a model for DNA melting, modelled as DNA unzipping, 
which provides a way to illustrate the role of the Boltz- 
mann factor in a venue well-known to the students. Sec- 
ond, we discussed the activity of helicase motor pro- 
teins in unzipping double-stranded DNA, for example, 
at the replication fork, which is an example of a Brown- 
ian ratchet. Finally, we treated force generation by actin 



polymerization, which is another Brownian ratchet, and 
for which we can determine how the velocity of actin 
polymerization depends on actin concentration and on 
load. In each of these examples, building on an ear- 
lier coverage of biased random walks, biology and pre- 
medical students in an introductory physics sequence at 
Yale were lead to the realization that a physics-based 
approach permits a deeper understanding of a familiar 
biological phenomenon. 
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